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ABSTRACT 

We report the detection of broad absorptions due to Si IV AA4088-4116 in the 
Luminous Blue Variable (LBV) AG Garinae during its last hot phase (2001-2003). 

Our NLTE spectral analysis, with the radiative transfer code GMFGEN, revealed 
the photospheric nature of these lines predicting, however, much narrower and 
deeper absorption profiles than observed. Using a recently-developed code to 
compute synthetic spectra in 2D geometry allowing for the effects of rotation, we 
could match the broad absorptions with a high projected rotational velocity of 190 
± 30 kms“^ on 2001 April. Analysis of spectra obtained on 2002 March and 2003 
January, when the star was cooling, yielded to a projected rotational velocity of 
110 ± 10 km s“^ and 85 ± 10 km s“^, respectively. The derived rotational velocities 
are proportional to as expected from angular momentum conservation. We 
discuss the effects of such high rotation on the spectral analysis of AG Gar, and 
on the wind terminal velocity. Our results show direct spectroscopic evidence, 
for the first time, that a LBV may rotate at a significant fraction of its break-up 
velocity. Thus, AG Gar (and possibly other LBVs) is indeed close to the LG 
limit, as predicted by theoretical studies of LBVs. 
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1. Introduction 

Massive stars have a tremendous impact on both the chemical and dynamical evolution 
of a galaxy. They input energy and momentum, and eject their nucleosynthesis products into 
the interstellar medium (Freyer et ah 2003). A detailed knowledge of massive star evolution 
is thus crucial to understanding their host galaxies, especially in the early Universe when 
massive stars are dominant (Matteucci & Calura 2005). 

Over the last decade, signihcant efforts have been made to include the effects of rotation 
consistently in the evolutionary models for massive stars (Meynet & Maeder 1997). It has 
been found that the effects of rotation are far from negligible, especially for the most mas¬ 
sive stars with high initial rotational velocity. As a consequence of rotation, the evolutionary 
tracks on the HR diagram, and the duration of the evolutionary stages (main sequence. Lu¬ 
minous Blue Variable - hereafter LBV, and Wolf-Rayet phases), are substantially modihed. 
For a recent review we refer to Maeder & Meynet (2000b). 

Rotation is also expected to play an important role in both LBV variability and LBV 
eruptions (Langer et ah 1999; Maeder & Meynet 2000a; Aerts et al. 2004). Even before 
reaching the Eddington limit, these stars may lose large amounts of material if they have a 
signihcant rotational velocity compared to their critical rotational velocity. Stellar rotation, 
during the LBV phase, can have an impact on the surrounding circumstellar medium, shaping 
their ejected nebulae with the characteristic bipolar geometry found around many of these 
stars (Humphreys & Davidson 1994; Maeder & Desjacques 2001; Dwarkadas & Owocki 2002). 

In this Letter we present high-resolution, high S/N optical spectra of the LBV AG 
Garinae. AG Gar has pronounced S Dor type variations, changing its effective temperature 
from 9,000 to 25,000 K (Leitherer et al. 1994; Stahl et al. 2001). AG Gar is considered the 
prototype of the LBV class and has been extensively observed, from UV to radio wavelengths, 
with photometric, spectroscopic and polarimetric techniques over several decades. It also 
has a bipolar nebula which is seen at various wavelengths (Nota et al. 1992; Voors et al. 


^Based on observations made at the 1.52 m ESO telescope at La Silla, Chile and at the 8 m VLT/ESO at 
Paranal, Chile 
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2000; Duncan & White 2002) that was ejected approximately 10"^ years ago (Lamers et ah 

2001 ). 

Here we report the detection of broad absorptions associated with the high-excitation 
lines of Si IV AA4088-4116 in AG Car during its last hot phase (2001-2003). These data are 
described in Sec. 2, while in Sec. 3 we present the results of our NLTE spectral modeling. 
Our modeling reveals the photospheric nature of the Si IV lines and allow us to derive, for 
the first time, the rotational velocity (urot) of AG Car and how it varies during the S Dor 
cycle. We show that Vrot is high compared to its critical rotational velocity {vcrit)- In Sec. 
4 we discuss some of the effects expected to be present in fast-rotating stars, such as the 
reduction in the wind terminal velocity and the presence of asymmetric winds. The high Vrot 
derived has important implications not only for AG Car and its wind, but for the general 
understanding of the LBV phase and the subsequent evolutionary stages of a massive star. 


2. Observations 

The spectroscopic observations of AG Car analyzed in this work were obtained as part 
of an observational program to follow the star over a complete S Doradus cycle with a 
dense time sampling and at high resolution. The resnlts provided by this campaign will be 
published in a forthcoming paper (Groh et al. in prep.). 

During 2001 and 2002, we collected AG Car spectra at the 1.52 m telescope of the Eu¬ 
ropean Southern Observatory (ESO) in La Silla, Chile, with the hber-fed, bench-mounted 
echelle spectrograph EEROS (Kaufer et al. 1999). To rednce the data we utilized an au¬ 
tomatic pipeline reduction system; further details are given in Stahl et al. (1999). The 
conhguration used on the EEROS spectrograph provided a resolntion of 48,000 over the re¬ 
gion 3600-9200 A, with increasing signal-to-noise ratio (S/N) from blue to red. We observed 
AG Car with EEROS on 2001 April 12 and 2002 March 17. 

In addition to the EEROS spectra, we retrieved spectra of AG Car from the ESO 
public-data archive^. They were obtained on 2003 Jannary 11 with the UVES spectrograph 
attached to the 8m VLT-Kneyen at ESO, under ESO Director Discretionary Time Program 
266.D-5655 (Bagnulo et al. 2003). The resolution provided by the UVES spectrograph was 
80,000 in the wavelength range 3000-10400 A, with an exceptional S/N of about 400. As 
with the EEROS spectra, we utilized the pipeline-processed data provided online. Table 1 
summarizes the high-resolution spectroscopic observations utilized in this Letter. 


^http://archive.eso.org 
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In Figure 1 we show (in velocity space) the absorption prohles of Si IV AA 4088-4116 
detected in AG Car during its last hot phase. From bottom to top, we show the 2001 April, 
2002 March and 2003 January spectra. The spectra were corrected to the heliocentric frame 
and the systemic velocity of lOkrns”^ (Stahl et ah 2001) was subtracted. The same broad 
Si IV absorption prohles can also be seen in high-resolution FEROS spectra that we have 
taken in 2001 January, 2001 June and 2002 July, that are omitted here due to their low S/N, 
and in lower resolution spectra. Hence, we conclude that the Si IV broad absorptions are 
present during the whole hot phase (i.e., they are not a transient feature). 


3. Modeling and results 

To analyze the observational spectra of AG Carinae we used CMFGEN (Hillier & Miller 
1998). CMFGEN assumes an outhow in steady-state, spherical symmetry and treats line 
and continuum formation in NLTE. Each model is specihed by the hydrostatic radius i?*, the 
luminosity L, the mass-loss rate M, the volume hlling-factor / and the chemical abundances 
Zi of the included species. The code does not solve the momentum equation of the wind, so 
a velocity law must be specihed. The velocity structure is described by a beta-law modihed 
at depth to match an hydrostatic structure at i?* (typically at a Rosseland depth of 100 in 
our models), so we also have to input the wind terminal velocity v^o and the steepness of the 
velocity law 13. We assumed a constant turbulent velocity of 20kms“^. A higher turbulent 
velocity would give a detectable redshift, which is not observed, of the emission lines with 
high opacity, such the stronger lines of H I and He I. 

While GMFGEN does not handle directly the effects of rotation in the radiative transfer, 
we accounted for its inhuence on the observed spectrum following the procedures described 
by Busche & Hillier (2005). For each GMFGEN model we computed the emerging spectrum 
in 2D geometry viewed from pole (0 degrees, which was equivalent to a non-rotating model), 
30 degrees and equator (90 degrees), for a given value of Vrot- As our goal is only to show 
that the broad absorptions are photospheric and broadened by rotation, we did not explore 
the effects of an equatorial or polar density enhancement on the line prohles. As shown by 
Busche & Hillier (2005), the hrst order impact will be on the relative line strengths, which 
would hrstly affect the derivation of the temperature and/or mass-loss rate. These effects 
will be investigated in a forthcoming paper. 

In Figure 2 we show the absorption prohles of Si IV A A 4088-4116 observed on 2001 
April compared with our best model. Based on the spectroscopic analysis, we obtained for 
this epoch L=1.0xl0® L©, R*=56 R©, T*=24,340 K, Voo=105kms“^, M=2.3xl0“^ M© yr“^, 
and f=0.15. The non-rotating model predicts signihcantly narrower and deeper absorption 
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profiles for both Si IV lines than the observations. As these Si IV lines are formed deep 
in the wind, the broadening effects due to the velocity held are negligible. We needed 
sini=190±30kms“^ to reproduce the observed lines’ shape and width, as can be seen 
analyzing the hts of model 269_2D^ and 269_2D_c. In spite of the noise of the data, we 
can clearly rule out a lower value of Vrot (e. g. a non-rotating model or model 269_2D_h with 
Vrot sin i =120 kms“^). 

To derive Vcru, and the surface escape velocity {vesc), we adopted the mass of AG Car 
as M = 35 Mq (Smith et ah 1994; Vink & de Koter 2002), which should be regarded as a 
lower limit (see discussion in Sec. 4 of Vink & de Koter 2002). Changes in M will affect the 
analysis, since Vcrit and Vesc are both proportional to and the Eddington parameter 

r depends on M~^. During 2001 April we calculated ncrit=220kms“^, following Maeder & 
Meynet (2000a), for a derived value of r=0.7. Assuming that we are seeing AG Car edge-on, 
the derived Vrot=190 kms“^ amounts to 0.86 Vcru- This value is exceptionally high and it is 
a lower limit, since higher values will be inferred for lower inclination angles i. For example, 
iii = 60°, we would have ^^ 0 ^= 220 kms“^ and Vrotlvcrit=^i he., the star would rotate nearly 
at Vcrit' 

The precise determination of Vrot/vcrit is beyond the scope of this Letter, since it relies 
on the accuracy of M, L, distance, reddening law, and i. However, for any reasonable change 
in one or more of these parameters, a signihcant value of Vrot/vcrit would still be reached. 
Hence, the rotationally-broadened photospheric absorptions reported in this Letter make it 
implausible that AG Car is seen pole-on or from lower inclination angles, as Vrot would be 
arbitrarily high. Moreover, the kinematics and dynamics of its bipolar nebula show extremely 
low projected velocities, in the range ±50kms“^ (Nota et al. 1992), which is consistent with 
a two-lobe polar ejection viewed nearly from the equator. 

It is important to point out that our models with rotation predict undetectable changes 
on the computed prohles of H I, He I lines, and other strong lines on AG Car’s spectrum. 
As they are formed far from the photosphere, the influence of rotation is very low, due to 
the r~^ effect (conservation of angular momentum). 

The spectra obtained during 2002 March have a S/N higher than the 2001 April dataset. 
In Figure 3 we present the best model for this epoch, compared with the observational 
data. The physical parameters of the star derived from the best hts of CMFGEN were 
L=l.OxlO®L 0 , R*=87R0, T*=19,650 K, Voo=195kms-^ M=4.7xl0-® M 0 yr-^ andf=0.25. 
The non-rotating model is inconsistent with the observed lines’ width and strength. To match 
the width and the shape of the lines we had to use a 2D model with sini=110kms“^. 
The high S/N allowed us to constrain the uncertainties to ±10kms“^. In spite of the superb 
hts to the line cores and red wings, we can also see in Figure 3 a high-velocity absorption in 
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the blue wings, which could not be reproduced by our models. Models with higher M could 
not reproduce this feature either. 

We present in Figure 4 the 2003 January observations around the Si IV lines com¬ 
pared with our best model. At this epoch AG Car was getting colder, as conhrmed by 
our derived stellar temperature of 17,940 K. Thanks to the high-resolution and high S/N 
of the VLT data, we could still clearly observe broader Si IV absorptions compared with 
the predictions of the non-rotating model, with parameters L=l.OxlO®L 0 , R*= 1 O 3 R 0 , 
T*=17,940 K, noo=150kms“^, M=6.0xl0“^ M 0 yr“^, andf=0.25. Running a 2D model with 
Urot sin i=85 km s“^, we could achieve a superb £t to the absorptions’ width and strength. A 
model with sini= 100 kms“^ has broader profiles than the observed ones, especially for 
Si IV A4116. Hence, we could constrain sini=85 ±10kms“^ in this epoch. 

In Table 2 we summarize the results obtained for the 3 epochs analyzed in this Letter. 
The derived luminosity agrees within 20% with previous works (Smith et ah 1994), while the 
mass-loss rate is clearly variable. A full discussion will be presented in a forthcoming paper. 
A key result is presented in the last column, where we show the derived value of Vrot R* ■ This 
physical parameter should be constant over time due to angular momentum conservation, 
as the stellar mass is practically invariant on this timescale. Analyzing the values derived in 
this work, we could verify that the rotational velocities obtained are consistent, taking the 
errors into account. This reinforces our finding that the broadening detected in the Si IV 
photospheric lines is a consequence of rotation. 


4. Discussion 

Rotation is believed to play a significant role dnring the evolution of massive stars, 
especially during the short-lived, transitional stages like the LBV phase, where the stellar 
and wind physical parameters change dramatically on short timescales. Rotation has also 
been invoked to explain the LBV phenomenon and the origin of asymmetric winds fonnd in 
some LBVs, like the snpermassive star rj Car (Smith et ah 2003; van Boekel et al. 2003), 
and related objects like B[e] stars (Zickgraf et al. 1985). For AG Car, observational evidence 
for wind asymmetry is provided by polarization data (Leitherer et al. 1994; Schulte-Ladbeck 
et al. 1994; Davies et al. 2005) and may have a significant impact on the spectroscopic 
analysis. Snch effects are cnrrently being explored by ns. 

Rotation, and its effects in the wind, have to be taken into account in any theory used 
to explain how the physical properties change during the S Doradus cycle for AG Car and 
other LBVs. One of the most noticeable effects of a high rotation, for example, would be a 
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significant decrease in the wind terminal velocity (Friend & Abbott 1986). Indeed, based on 
our derived physical parameters, we can calculate that during the hottest epoch (April 2001) 
nesc=290 kms“^, which includes the reduction due to the radiation pressure. The derived 
wind terminal velocity of 105kms“^ is incompatible with that surface escape velocity and 
the predicted value of Voo/vesc ~ 1.3 —2.6 for hot supergiants (Tamers et ah 1995). However, 
if we account for the reduction effect that a high rotational velocity has on the effective 
gravity geff (Maeder & Meynet 2000a), and hence on the surface escape velocity (dehned as 
Vesc = {‘^9eff we would achieve Vesc ~ 40 — 70kms“^ and Voo/vesc ~ 1-5 — 2.6. 

The high rotational velocity also explains why was lower when the star was hotter 
(2001 April) compared to when the star was cooling (2002 March and 2003 January). At 
a hrst glance, this result disagrees with the predicted temperature dependence of Voo for 
line-driven winds (Vink et ah 2000). However, the effect of rotation on reducing the wind 
terminal velocity is higher during the hot phase, which can be translated to higher values of 
Vrot/vcritj as showu in Table 2. This can be understood in terms of the R~^ dependence on 
the rotational velocity, while the critical rotational velocity depends on 

So far, a high rotational velocity has never been observed directly (i. e. in the stellar 
spectrum) in any active LBV, in spite of the intense theoretical studies made in recent years. 
Our spectroscopic analysis shows that AG Car, a well-known prototype of the LBV class, 
rotates at a signihcant fraction of its break-up velocity. Interestingly, the behavior of the 
rotational velocity as a function of the stellar radius during 2001-2003 is consistent with 
a R~^ dependence (see Sec. 3), as expected from angular momentum conservation. Our 
modeling also conhrms that AG Car is very luminous (L=1.0xl0®L©), and therefore it lies 
close to the so-called TO limit (Maeder & Meynet 2000a). Thus, it is likely that rotation, 
coupled with a high luminosity, plays a key role during the S Dor cycles, the LBV giant 
eruptions and the subsequent late evolutionary stages of a massive star. 
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Hillier and the University of Pittsburgh for hospitality and support for this work. We thank 
an anonymous referee for useful suggestions. 
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Table 1: Summary of AG Car spectroscopic observations 


Epoch 

Telescope 

Spectral Range (A) 

R 

2001 Apr 12 

1.52m ESO 

3600-9200 


2002 Mar 17 

1.52m ESO 

3600-9200 


2003 Jan 11 

8m ESO (VLT) 

3000-10400 
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Table 2: Derived physical parameters of AG Car during 2001-2003, assuming i = 90° 


Epoch 

T* 

M (lO-^Moyr-i) 

Voo (kms“^) 

Nrot (kms“^) 

cvit 

R*(R0) 

v^otR* (lO'^Rokms”^ 

2001 Apr 12 

24,340 

2.3 

105 

190 ± 30 

0.86 

56 

10.7 ± 2.2 

2002 Mar 17 

19,650 

4.7 

195 

110 ± 10 

0.63 

87 

9.5 ± 1.7 

2003 Jan 11 

17,940 

6.0 

150 

85 ± 10 

0.57 

103 

9.3 ± 1.8 


We assume a conservative 10% error on R+, accounting for uncertainties on the distance, reddening law and 
luminosity. 
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Fig. 1.— High resolution, continuum-normalized spectra of AG Car, centered on Si IV 
4088.86 A (left panel) and on Si IV 4116.10 A (right panel). The spectra were offset for 
the sake of clarity. On both panels, we show from top to bottom data from 2003 January 
(-1-0.05 offset), 2002 March (no offset) and 2001 April (-0.15 offset). Signihcant changes in 
the strength of the Si IV lines, related to changes in the effective temperature of the star, 
can be seen. 
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Fig. 2.— High-resolution observed absorption profiles (full line) at Si IV 4088 (left panel) 
and Si IV 4116 (right panel). A set of CMFGEN models, computed without rotation (dot- 
dashed), and with a projected rotational velocity of 120kms“^ (dotted line), 180kms“^ 
(double dot-dashed) and 220kms“^ (dashed line). For the sake of clarity, at the right panel 
we only show the non-rotating, 180kms“^ and 220kms“^ models. In spite of the errors due 
to the low S/N, we can clearly rule out a rotational velocity lower than ~150kms“^ for AG 
Gar on this epoch. 
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Fig. 4.— Same as Fig. 2, but for the spectra gathered on 2003 January (full line). The CMF- 
GEN models were computed with (dashed) and without (dot-dashed) a projected rotational 
velocity of 85kms“^. To estimate the uncertainties, the dotted line shows the CMFGEN 
model with a projected rotational velocity of lOOkms”^. 













